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One limit to the fidelity of quantum logic operations on trapped ions arises from heating of the 
ions' collective modes of motion. Sympathetic cooling of the ions during the logic operations may 
eliminate this source of errors. We discuss benefits and drawbacks of this proposal, and describe 
possible experimental implementations. We also present an overview of trapped-ion dynamics in 
this scheme. 



I. INTRODUCTION 



One of the most attractive physical systems for generating large entangled states and realizing a quantum computer 
101 is a collection of cold trapped atomic ions ||^. The ion trap quantum computer stores one or more quantum 
bits (qubits) in the internal states of each trapped ion, and quantum logic gates (implemented by interactions with 
externally applied laser beams) can couple qubits through a collective quantized mode of motion of the ion Coulomb 
crystal. Loss of coherence of the internal states of trapped ions is negligible under proper conditions but heating of 
the motion of the ion crystal may ultimately limit the fidelity of logic gates of this type. In fact, such heating is 
currently a limiting factor in the NIST ion-trap quantum logic experiments pfl. 



Electric fields from the environment readily couple to the motion of the ions, heating the ion crystal If the 

ion trap is much larger than the ion crystal size, we expect these electric fields to be nearly uniform across the crystal. 
Uniform fields will heat only modes that involve center-of-mass motion (COM motion), in which the crystal moves as 
a rigid body. Motional modes orthogonal to the COM motion, for instance the collective breathing mode, require field 
gradients to excite their motion. The heating of these modes is therefore suppressed [Q. However, even if quantum 
logic operations use such a "cold" mode, the heating of the COM motion can still indirectly limit the fidelity of logic 
operations. Since the laser coupling of an internal qubit and a motional mode depends on the total wavepacket spread 
of the ion containing the qubit, the thermal COM motion can reduce the logic fidelity ||,^. 



In this paper, we examine sympathetic cooling |g] in a particular scheme for which we can continuously laser-cool 
the COM motion while leaving undisturbed the coherences of both the internal qubits and the mode used for quantum 
logic. In this method, one applies continuous laser cooling to only the center ion of a Coulomb-coupled string of an odd 
number of ions. One can address the center ion alone if the center ion is of a different ion species than that composing 
the rest of the string |^ . Alternatively, one can simply focus the cooling beams so that they affect only the center ion. 
In either case, the cooling affects only the internal states of the center ion, leaving all other internal coherences intact. 
If the logic operations use a mode in which the center ion remains at rest, the motional coherences in that mode are 
also unaffected by the cooling. On the other hand, the sympathetic cooling keeps the COM motion cold, reducing 
the thermal wavepacket spread of the ions. In the following, we will discuss the dynamics of an ion string in which all 
ions are identical except the center ion, assuming heating by a uniform electric field. Our results give guidelines for 
implementing the sympathetic cooling scheme. Similar results would apply to two- and three-dimensional ion crystals 
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II. AXIAL MODES OF MOTION 



We consider a crystal of N ions, all of charge q, in a linear RF trap [ |lOyil[| . The linear RF trap is essentially an RF 
quadrupole mass filter with a static confining potential along the filter axis z. If the radial confinement is sufhcicntly 
strong compared to the axial confinement, the ions will line up along the z-axis in a string configuration ||l0| , [rf[| . There 
is no RF electric field along z, so we can write the axial confining potential as (f>{z) — for ao a constant. The 

potential energy of the string is then given by 

f ^ q2 ^ 1 

V{z, z„) = -qao E + ^ E ^-3^ (1) 

for Zi the position of the ith ion in the string (counting from the end of the string). The first term in the potential 
energy expresses the influence of the static confining potential along the z-axis, while the second arises f rom the 
mutual Coulomb repulsion of the ions. For a single ion of mass m, the trap frequency along z is just uj^ = \J qa^^/m. 

We can compute the equilibrium positions of the ions in the string by minimizing the potential energy of Eq. 
Defining a length scale ihy (.^ — q/(47reoao) and normalizing the ion positions by Ui = Zijl gives a set of equations 
for the Ui as 
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= 0, i = f .. .TV (2) 



which has analytic solutions only up to = 3. Steane [Q and James ||l4| have computed the equilibrium positions 
of ions in strings with TV up to 10. The potential energy is independent of the mass, so the equilibrium positions of 
ions in a string are independent of the elemental composition of the string if all the ions have the same charge. 



In a real ion trap the ions will have some nonzero temperature and will move about their equilibrium positions. If 
the ions are sufficiently cold, we can write their positions as a function of time as Zi(i)) — £ui + qi(t), where qi{t) is 
small enough to allow linearizing all forces. We specialize to the case of an odd number of ions N , where all ions have 
mass TO, except for the one at the center of the string which has mass M . The ions arc numbered 1 . . . A^, with the 
center ion labeled by Uc = {N + l)/2. Following James flj], the Lagrangian for the resulting small oscillations is 
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where 
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We define a normalized time as T = LOzt- In treating the case of two ion species, we write ji — M/m for the mass ratio 
of the two species and normalize the amplitude of the ion vibrations qi{t) as Qi = ft^qao, i 7^ n-c, Qn^ — <ln^^/^aQil. 
The Lagrangian becomes 
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where 



Aij i, j ^ Uc 
A'i, = I Aij/^ i or j ^ He, j 



(7) 



generalizing the result of James |lj 



The Lagrangian is now cast in the canonical form for small oscillations in the coordinates Qi (t) . To find the normal 
modes, we solve the eigenvalue equation 



A' • = CfcV'=) k=l...N 



(8) 



for the frequencies Cfc and (orthonormal) eigenvectors v'^'^-' of the iV normal modes. Because of our normalization of 
the Lagrangian (|^), the Cfc are normalized to ujz and the v^''"' are expressed in terms of the normalized coordinates 
Qi{t). In terms of the physical time the frequency of the fcth mode is Ck^z- If the fcth mode is excited with an 
amplitude C, we have 



(9) 
(10) 



in terms of the physical coordinates qi{t). 



We can solve for the normal modes analytically for A'^ = 3. Exact expressions for the normal- mode frequencies are 
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normalized to w^. The mode eigenvectors are 



(21 + v/441 - 34/i + 169/i2) 



v^^^ (1,0,-1) 
.-(3)=A.3(^l,N^(13-5C|),l) 



(11) 
(12) 

(13) 



(14) 
(15) 
(16) 



in terms of Qi{t). Here Ni,N2, and N3 are normalization factors. In the case of three identical ions (/i — 1), we can ex- 
press the mode eigenvectors in terms of the Qi(i) as w^^) = (1, 1, 1)/V3, z/^^) = (1, 0, — 1)/'\/2, and w^^^ = (1, — 2, 1)/V6. 
The mode eigenvectors, in this special case, also give the ion oscillation amplitudes in terms of the physical coordinates 
qi{t). For three identical ions, then, pure axial COM motion constitutes a normal mode. (This result holds for an 
arbitrary number of identical ions.) We also note that the center ion does not move in mode #2; hence the frequency 
and eigenvector of mode #2 are independent of /i. For any odd number N of ions there are {N — l)/2 modes for 
which the center ion does not move. These modes will likewise have frequencies and eigenvectors independent of /z. 



Moreover, they have v. 



(fc) 



,(fe) 



and so they are orthogonal to the COM motion and do not couple to uniform 



electric fields. The center ion moves in the other (A^ -I- 1)/2 modes, and unless n = I, each of these {N + l)/2 modes 
has a component of axial COM motion and therefore couples to uniform electric fields. 



For N = 5 and higher, the normal mode frequencies depend on in a complicated way. However, it is easy to 
find the frequencies numerically. Fig. 1 shows the mode frequencies for TV = 3, 5, 7, and 9 as a function of fi for 
0.01 < fi < 100. The modes are numbered in order of increasing frequency (at = 1), and are normalized to tOz- 
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In each case, the lowest-lying mode has all ions moving in the same direction and consists of pure COM motion for 
jj, = I. The even- numbered modes correspond to the {N — l)/2 modes for which the center ion does not move. Their 
frequencies are therefore independent of fi. For both very large and very small fi the modes pair up, as shown in Fig. 
1. For each pair there is some value /i > 1 for which the modes become degenerate. The relative spacing between 
modes in a pair is also smaller in the large-/x limit than in the small-/x limit. 

In selecting a normal mode of motion for logic operations, we want to ensure that the mode is well resolved from 
all other normal modes. However, when modes are nearly degenerate, as for fi <^ 1 and /i ^ 1, transfer of energy 
can occur between the modes in the presence of an appropriate coupling, for instance if the static confining potential 
contains small terms of order This coupling can lead to a loss of coherence of the logic mode. Also, the need 

to resolve the logic mode from a nearby spectator mode can force a reduction in gate speed. These effects limit the 
usefulness of the sympathetic cooling scheme for /z very large. Evidently it is best to use a cooling ion that is of the 
same mass or lighter than the logic ions. In this case mode ^2 is well-separated from all other modes, as shown in 



We now consider the motion of the ions transverse to the z-axis. The ions experience an RF potential 
XCOs{ftt){x'^ — j/^)/2 for a suitable choice of axes x and y perpendicular to z, where il is the frequency of the 
RF field and x is a constant. The static confining potential can be written (qao/2)(z^ — ax'^ — (1 — a)y'^) at the 
position of the ions (with a a constant), so there is also a transverse static electric field. To analyze the ion motion, 
we work in the pseudopotential approximation, in which one time-averages the motion over a period of the RF drive 
to find the ponderomotive force on the ion. If the static potential is negligible, the RF drive gives rise to an effective 
transverse confining potential of ^muj^Q{x^ + y^), where w^o — ^x/iV^^m) for an ion of mass m. If we include 
the effects of the static field, the transverse potential becomes ^mluj-^x'^ + i^Jyy^), where iVx — i^roy/^ — Oi<^z/^rO^ 
ujy = WrO\/l ~ (1 ~ oi)ljjHu)1q. Below we will assume a — 1/2, so that ujy — ujx- In any case, the transverse potential 
is that of a simple harmonic oscillator, as we saw also for the axial potential. However, the transverse potential 
depends directly on the ion's mass, so the center ion of a string feels a different trap potential than the others for 



We define e = t^ro/^^z, so that lu^ — LUzy^e^ — 1/2. Then the normalized Lagrangian for the motion along x is 



Fig. 1. 



III. TRANSVERSE MODES OF MOTION 




(17) 



where Xi — Xi^^qao for i ^ ric and X, 



Xi^/qooji are normalized ion vibration amplitudes along x. Here 




(18) 



and 





(19) 
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We can describe the normal mode frequencies and oscillation amplitudes in terms of the eigenvectors and eigenvalues 
of B'^j^ just as for the axial case above. The normalizations of the time and position coordinates remain the same as 
in the axial case. 



In the previous section, we assumed that the radial confinement of the ions was strong enough that the configura- 
tion of ions in a string along the z-axis was always stable. However, for sufficiently small e, the string configuration 
becomes unstable. The stable configurations for different values of e can be calculated and several of these 

configurations have been observed for small numbers of ions Rather than review the theory of these configu- 

rations, we will simply find the range of validity of our small-oscillation Lagrangian for the string configuration. The 
string will remain stable for all e greater than some = es(/Lt); also varies with N. On the boundary between 
stable and unstable regions, the frequency of some mode goes to zero. Recalling that the determinant of a matrix 
is equal to the product of its eigenvalues, we see that es(/i) is the maximum value of e satisfying det-B'(e, /i) = 
for /i fixed. Fig. 2 shows es(/i) as a function of fi for 3, 5, 7, and 9 ions. In each case, there is a cusp in es(/i) 
corresponding to the crossing of the two largest solutions to det B'{e, fi) = 0. The position of the cusp varies with the 
number of ions, but lies between /i = 0.1 and /i = 1 for < 9. Only the cusp for = 3 is clearly visible in Fig. 
2, but numerical study indicates the presence of a cusp for all four values of N. For fi greater than the value at the 
cusp, e < esifJ-) corresponds to instability of the zigzag mode, so that the string breaks into a configuration in which 
each ion is displaced in the opposite direction to its neighbors p^ , |l7| . For /i smaller than the value at the cusp, is 
independent of /i. In this regime, e < creates an instability in a mode similar to the zigzag mode, except that the 
center ion remains fixed. 



We can proceed to calculate the frequencies of the transverse modes for values e > esifJ.). Again, these frequencies 
are normalized to the axial frequency of a single ion of mass m. Fig. 3 shows the transverse mode frequencies for 3, 
5, 7, and 9 ions as a hmction of /i, where e is taken equal to l.les(/i). The modes are numbered in order of increasing 
frequency at fj, — 1 (all ions identical). In this numbering scheme, the central ion moves in odd- numbered modes 
but not in even-numbered modes. The frequencies of the even-numbered modes appear to depend on /i because 
they are calculated at a multiple of Ssip), for constant e these frequencies are independent of fi. The cusps in the 
mode frequencies in Fig. 3 arise from the cusps of es{fJ.) at the crossover points between the two relevant solutions 
of det B' — 0. Mode frequencies plotted for a constant value of e do not exhibit these cusps. As in the case of axial 
motion, the mode frequencies form pairs of one even- and one odd-numbered mode for small fj,. However, for large fj, 
all but one of the transverse modes become degenerate. The only nondegenerate transverse mode in this case is the 
zigzag mode. In general, the modes are most easily resolved from their neighbors for /i = 1, as in the case of axial 
motion. Increasing e reduces the frequency spacing between nearly degenerate modes. At e = l.les{fJ,) and fJ, = 1, for 
instance, the fractional spacing between the cold transverse mode of 3 ions and its nearest neighbor is 0.20, but for 
e = 1.5es(/i) the same spacing is 0.09. 

The near-degeneracy of the modes for large or small /i and for e/eg significantly greater than 1 limits the usefulness 
of these modes because of possible mode cross-coupling, just as for the axial modes. Resolving a particular transverse 
mode requires operating the trap near the point at which the string configuration becomes unstable, i.e., e near eo(/u). 
In this regime, the collective motion of the ions is quite sensitive to uncontrolled perturbations, which may pose 
significant technical problems for using a transverse mode in quantum logic operations. 



IV. MODE HEATING 

Stochastic electric fields present on the ion trap electrodes, for instance from fluctuating surface potentials, can 
heat the various normal modes of motion incoherently. For ion trap characteristic dimension dtrap much larger than 
the size of the ion crystal dions, these fields are approximately uniform across the ion crystal, so they couple only 
to the COM motion. The (iV — l)/2 even-numbered modes are orthogonal to the COM motion, so they are only 
heated by fluctuating electric field gradients. The heating rates of these modes are reduced by a factor of at least 
{dions / dtrapY ^ 1 as Compared to the heating of the other modes 0. In the following, therefore, we will neglect the 
effects of fluctuating fleld gradients, so that the even-numbered modes do not heat at all. 
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The analysis of sections 2 and 3 shows that the motion of a crystal of N ions is separable into the 3N normal 
modes, each of which is equivalent to a simple harmonic oscillator. Hence we can quantize the crystal motion by 
quantizing the normal modes. The fcth normal mode gives rise to a ladder of energy levels spaced by hC,k^z, with 3A^ 
such ladders in all. If we now write the uniform electric field power spectral density as Se{ijj), we can generalize the 
result of ||l8| to give 



nk = — — TT 



(k) 



N 



(fc) 



(20) 



for the heating rate of the fcth mode, expressed in terms of the average number of quanta gained per second. Recall 

(k) 

that Uj is the oscillation amplitude of the ith ion in the fcth normal mode, expressed in the normalized coordinates. 
It is useful to normalize the heating rate in equation ( pO| ) to the heating rate of the lowest-lying axial mode of a string 
of identical ions. This normal mode consists entirely of COM motion and we write v'^^^'^ — for all ions. The 

normalized heating rate of the fcth mode is then 



rifc 



ncoM 



(fc) N 



(21) 



where we have assumed that the spectral density S'£;(w) is constant over the frequency range of the normal modes, 

i.e., SEii^z) = SEiCkUJz)- 



Fig. 4 shows plots of the normalized heating rates of the axial modes for = 3, 5, 7, and 9 as a function of ^. Fig. 
5 is the same, but for the transverse modes, with e = Lle^. The numbering of modes on the plots of heating rate 
matches the numbering on the corresponding plots of mode frequency (Figs. 1 and 2). 



In both axial-mode and transverse-mode plots, the even-numbered modes have the center ion at rest, while the 
center ion moves for all odd-numbered modes. We see from Figs. 4 and 5 that the modes for which the center ion is 
fixed can never heat, while all the other modes always heat to some extent for /i 7^ 1. We will refer to these modes 
as "cold" and "hot" modes, respectively. If the ions are identical, only the modes with all ions moving with the 
same amplitude (COM modes) can heat. There are three such modes, one along i, one along y, and one along z. In 
interpreting Figs. 4 and 5, it is important to recall that the normalized heating rate defined in Eq. ( pT| ) is inversely 
proportional to the mode frequency. For instance, the ^-dependence of the heating rate of the highest-frequency 
transverse mode can be largely ascribed to variations in the mode frequency, rather than to changes in the coupling 
of the mode to the electric field. 



V. PROSPECTS FOR SYMPATHETIC COOLING 



Heating reduces logic gate fidelity in two ways. The logic mode itself can be heated, but by choosing a cold mode, we 
can render this effect negligible. On the other hand, the Rabi frequency of the transition between logic-mode motional 
states depends on the total wavepacket spread of the ion involved in the transition . Heating on modes other 
than the logic mode can thus lead to unknown, uncontrolled changes in this Rabi frequency, resulting in overdriving 
or underdriving of the transition. The purpose of sympathetic cooling is to remove this effect by cooling the center 
ion and thus all hot modes. 



For sympathetic cooling to be useful, we must find a cold mode suitable for use in quantum logic. The cold mode 
must be spectrally well separated from any other modes in order to prevent unwanted mode cross-coupling. We can 
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use the lowest-lying cold axial mode as the logic mode for < 3. In this mode, called the breathing mode, the center 
ion remains fixed and the spacings between ions expand and contract in unison. Unless the trap is operated very close 
to the instability point of the string configuration, the breathing mode is better separated from its neighbors than are 
any of the cold transverse modes. For ^ 3 any cold mode, cither axial or transverse, is nearly degenerate with a hot 
mode. In this regime one must make a specific calculation of mode frequencies in order to find the best-resolved cold 
mode. Even so, the cold axial modes are again better separated from their neighbors than are the cold transverse 
modes, except for e very close to es(/i). It seems best to select a cold axial mode as the logic mode in most cases. 

By selecting our laser-beam geometry appropriately, we can ensure that the Rabi frequency of the motional tran- 
sition on the axial mode used for logic depends chiefly on the spread of the ion wavepacket along z. In this case, 
heating of the axial modes will affect logic-gate fidelity, but heating of the transverse modes will have little effect. 
If the mass of the central ion is nearly the same as that of the others (/x « 1), only the lowest axial mode will heat 
significantly, and we can continuously cool this mode by cooling only the central ion, ensuring that all ions remain in 
the Lamb-Dicke limit [||. If [i is not near 1, we must cool all (N l)/2 hot modes (again by addressing the central 
ion) to keep all ions in the Lamb-Dicke limit. 

The analysis above indicates that, all other things being equal, we are best off if our substituted ion is identical 
to, or is an isotope of, the logic ions. However, sympathetic cooling can still be useful if the two ion species have 
different masses. For example, we can consider sympathetic cooling using the species ^Be+ and ^''Mg+. Linear traps 
constructed at NIST have demonstrated axial secular frequencies of over 10 MHz for single trapped ^Be"'" ions. For 
three ions with ^^Mg+ as the central ion, uj^ (Be+) = 27r X 10 MHz yields a spacing of 1.6 MHz between the cold 
axial breathing mode and its nearest neighbor. If we reverse the roles of the ions (u;^(Mg+) = 27r X 10 MHz), the 
spacing increases to 6.2 MHz. The transverse modes are much harder to resolve from each other. For three ions 
with ^^Mg+ in the center, we require WroCBe"*") = 27r x 27.6 MHz to obtain e = l.lcs, and the spacing between 
the cold transverse zigzag mode and its nearest neighbor is only 560 kHz. Reversing the roles of the ions, we find 
e — l.les at tj^o (Mg+) = 27r x 14.7 MHz with a spacing of 1.1 MHz. For this combination of ion species, the cold 
axial breathing mode seems most appropriate for logic. For a string of 3 or 5 ions, sympathetic cooling would require 
driving transitions on 2 or 3 axial- mode sidebands, respectively. From this example we see that sympathetic cooling 
can be useful even for ion mass ratios of nearly 3 to 1. 



VI. CONCLUSION 

We have investigated a particular sympathetic cooling scheme for the case of an ion string confined in a linear RF 
trap. We have numerically calculated the mode frequencies of the axial and transverse modes as functions of the 
mass ratio and trap anisotropy e for 3, 5, 7, and 9 ions. We have also calculated the heating rates of these modes 
relative to the heating rate of a single ion, assuming that the heating is driven by a uniform stochastic electric field. 
The results indicate that the scheme is feasible for many choices of ion species if we use a cold axial mode as the 
logic mode. The optimal implementation of the scheme employs two ion species of nearly equal mass. However, a 
demonstration of sympathetic cooling using ^Be''" and ^^Mg+ appears well within the reach of current experimental 
technique. 
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Fig. 1 




FIG. 1. Normalized axial mode frequencies as a function of n for (a) 3, (b) 5, (c) 7, and (d) 9 ions. 
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Fig. 2 
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FIG. 2. Trap anisotropy at instability of the string configuration as a function of /i for 3, 5, 7, and 9 ions. Arrows indicate 
the cusps discussed in the text. 
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Fig. 3 




FIG. 3. Normalized frequencies of the transverse modes as a function of with e = l.leo(/u) for (a) 3, (b) 5, (c) 7, and (d) 9 
ions. 
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Fig. 4 




Fig. 5 




FIG. 5. Normalized heating rates for the transverse modes as a function of n with e = l.leo(ju) for (a) 3, (b) 5, (c) 7, and 
(d) 9 ions. 
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